Decoherence of Rabi oscillations in a single quantum dot 
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We develop a realistic model of Rabi oscillations in a quantum-dot photodiode. Based in a multi- 
exciton density matrix formulation we show that for short pulses the two-level models fails and 
higher levels should be taken into account. This affects some of the experimental conclusions, such 
as the inferred efficiency of the state rotation (population inversion) and the deduced value of the 
dipole interaction. We also show that the damping observed cannot be explained using constant 
rates with fixed pulse duration. We demonstrate that the damping observed is in fact induced by 
an off-resonant excitation to or from the continuum of wetting layer states. Our model describes 
the nonlinear decoherence behavior observed in recent experiments. 
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Coherent manipulation of a quantum state is one of 
the tasks required in quantum computation and infor- 
mation processing. Semiconductor quantum dots (QDs) 
have been shown to be excellent candidates for the phys- 
ical implementation of these objectives. Rabi oscillations 
- temporal coherent oscillations of the population and its 
inversion in a two-level system driven by a strong reso- 
nant field - have been successfully demonstrated by dif- 
ferent groups using excitons in single QDs with different 
clever probing techniques fiiS42iiiS4S4L& as well as utilizing 
the biexciton population^ Zrenner et alA have developed 
a single self-assembled QD photodiode in which the pop- 
ulation inversion is probed by the photocurrent signal 
induced by a strong and carefully tuned optical pulse. 
In their device the pulse generates an electron-hole pair 
in the QD and with the help of an external gate volt- 
age, the electron and hole tunnel out into nearby con- 
tacts. This process generates a photocurrent signal that 
is a weakly disturbing probe of the coherent state of the 
system. Their results show Rabi oscillations that are 
damped with increasing area of the pulse for a fixed pulse 
duration. However, the mechanisms that produce this 
decoherence were unclear. 

We present a study of the dynamics of a single self- 
assembled quantum dot photodiode in the presence of an 
optical pulse. The electron and hole tunneling processes 
are introduced via a microscopic model of the structure 
which also includes the electron-hole interaction. The 
dynamics is described using a density matrix approach 
that incorporates dipole coupling to multi-exciton states 
and off-resonant excitation to states in the wetting layer 
(WL). Our model shows that for short pulses (of the or- 
der of a few ps, as used in experiments^) the two level 
system approximation fails for a 7r-pulse (the pulse neces- 
sary to invert the exciton population) and the biexciton 
state plays an important role in the dynamics. This re- 
sults in a frequency shift of the measured photocurrent 
oscillations, which can significantly affect the experimen- 
tally deduced value of the transition dipole moment. The 
efficiency estimate of the state rotation (population in- 
version) for a TT-pulse is also affected as the biexciton 
contributes to the photo-signal. We demonstrate that 
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FIG. 1: Schematic band structure and level configuration of 
QD photodiode. (a) An electromagnetic pulse creates an exci- 
ton in the dot, and an applied gate voltage forces the electron 
and hole to tunnel out, generating the photocurrent signal, 
(b) Schematic representation of processes and levels involved 
in this system: \x) and \b) are exciton and biexciton states. 



longer pulses minimize the biexciton contribution. Most 
importantly, we find that inclusion of excitations to WL 
states is essential to understand the decoherence observed 
in experiments. The damping due to coupling with the 
WL not only explains the observed shape of Rabi oscil- 
lations, but also gives a description of the background 
signal observed in Rabi-photodiodes.^° Although we fo- 
cus on the QD photodiodcf^ our model is also relevant for 
experiments on Rabi oscillations with optical readout.^ 

In Fig. n show the system and level configuration 
taken into account by our model Hamiltonian which can 
be written asii 
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Here, — Sx — fnu accounts for the detuning of the 
exciton with the laser energy huj, St — st — '2huj is the 
two photon biexciton detuning, ilx{t) = (0|/i • E{t)\x)/h, 
^b{t) — {x\fl-E{t)\b)/h, where the electric dipole moment 
fl describes the coupling of the excitonic transition to the 
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FIG. 2: Contour plot of the biexciton average occupation as 
function of pulse duration tp and pulse area Q. Notice long 
pulses do not produce biexciton population. Inset compares 
exciton and biexciton contribution for a 1 ps pulse. 

radiation field, and E{t) is the pulse amplitude which we 
assume to have a Gaussian shape, with full width at half- 
maximum (FWHM) tp (pulse length). 

For a strong pulse resonant with the exciton energy 
{5x = 0), the biexciton binding energy A = 3 meV is 
relatively small and, as can be seen in Fig. 13 cannot 
be neglected in the dynamical description of the system. 
There we show the average occupation of the biexciton 
state as a function of the pulse length tp and pulse area 
Q — U{t)dt, assuming that the biexciton has the 
same dipole moment of the exciton transitioniS [^b{t) — 
ilx{t) = ^{t)] and not including decay or dephasing. We 
show below that even after inclusion of these effects the 
biexciton population is important for short pulses. 

The dynamics of the reduced system is computed using 
a master equation in the Lindblad form.« 

|^4[i/o,p]+L(p), (2) 

where the first term on the right yields the unitary evo- 
lution of the quantum system and L(p) is the dissipative 
part of the evolution assuming the Markovian approxi- 
mation. We use Fj to describe all population decay rates 
of the level i, which we assume to have two different 
sources, one due to the spontaneous decay FJ'"^ given by 
the recombination rate (for the exciton state the spon- 
taneous decay time is known to be Tioc = 1/^^'^'^ ~ 1 ns 
for this kind of QD) and other F'^" which results when 
the particles (making the exciton or biexciton) leave the 
system by tunneling, so that F^ = F*"" -I- F™'^. 

We can estimate the tunneling rate F*"'^ for a 
single particle s (electron or hole) using the tun- 
neling Hamiltonian»i^'i^ The probability of OD to 
3D tunneling for a single particle isi^ F*^" — 
¥ Ea l(*.|^^Dot|*a)|' 5{E, - E^) , where and E, is 
the single particle wave function (for electron or hole) 



localized in the QD and its respective energy, while 
and Ea are correspondingly in the contacts. In 
we use plane waves for the in-plane motion, and 
the z-component was calculated using the linear poten- 
tial from the applied gate voltage and an exponential 
decay function in the barrier region^ ^'s is strongly 
localized in the quantum dot, and it is modelled as 
^'s(r, z) = -^^7=1- exp(— ^)x(z), the product of a ground 
state harmonic-oscillator function which describes the 
lateral motion in a quantum dot, and x(^)i the wave 
function of a square well potential. We obtain 

where m* is the effective mass of the particle s in the dot, 
L and Vs are the width and depth of the corresponding 
square well, and F is the electric field provided by the 
gate voltage. This is similar to the WKB tunneling but 
with well-characterized energy dependent prefactor. 

One can consider as a good first approximation that 
the electron-hole interaction produces a shift (the exciton 
binding energy) in the single particle energy level Es in 
the quantum dot, and use Eq. Oto obtain the rates with 
and without electron-hole interaction. Note that this es- 
timate of the carrier tunneling time r*"" = 1 /F^^" cannot 
be directly compared with the experimental linewidth of 
the photocurrent peak since that refiects the time for both 
particles to tunnel out of the system (even if dominated 
by the faster rate). Using parameters from Ref. ^3 and 
the assumption that Ve — 3Vk (where Ve (Vh) is the elec- 
tron (hole) depth of the square well profile), we obtain 
that the tunneling times for electron or hole are simi- 
lar when there is electron-hole interaction (this of course 
depends on QD parameters). After one of the particle 
tunnels, the remaining particle tunnels out faster, as it 
no longer experiences the electron-hole interaction. The 
rates we consider here, F|j."" for the exciton and Fj;"" for 
the biexciton, are the rates for all particles (electrons and 
holes) to leave the system, and they are obtained by solv- 
ing a separate density matrix equation including the dif- 
ferent tunneling rates for individual single-particle states. 
For the exciton state for instance, we solve a four dimen- 
sional density matrix which includes the vacuum, exciton, 
electron (when the hole leaves the system first) and hole 
states (when the electron leaves the system first). Solv- 
ing this density matrix one can evaluate the equivalent 
rate for a two level system that reproduces the multi- 
path process. We obtain F*""^ = where r*""^ ~ 12 
ps is the time for both particles to leave the system by 
tunneling, and is consistent with the photocurrent signal 
linewidth (~ 10 ps). We can consider this tunneling as a 
single process because the charged exciton, which could 
be created from a single particle state (after one particle 
leaves the QD), is ~ 4.6 meV out of resonance from the 
exciton state excitation for this kind of QD,^i and can- 
not be created during the pulse duration which is much 
faster than the single particle tunneling time. A similar 
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description applies to the biexciton state. This allows us 
to write a simple equation for the photocurrent signal in 
terms of the exciton and biexciton contributions as we 
now describe. 

As the photocurrent is a signal induced by the particles 
that tunnel out, when the next pulse arrives the system 
is already in the vacuum state |0) (meaning that they are 
different processes that account for a statistical average) , 
we can write the expression for the photocurrent as 
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Exciton with correction 
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where / is the repetition frequency of the pulse sequence 
(we use / = 82 MHz as in Zrenner's experiment*) and q 
is the electronic charge. One can also write this as 
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where r is the time at end of the pulse. Notice that the 
result is not a simple summation of exciton and biexci- 
ton contributions, but it is rather a mixture. This is ex- 
pected, as the biexciton occupation after the pulse even- 
tually decays to the exciton state, and contributes to the 
exciton part of the photocurrent. It is important to point 
out that even if we ignore the biexciton contribution to 
Eq. El it still exhibits a term not considered before^ the 
contribution during the pulse given by the second inte- 
gral in Eq. |S| This contribution can indeed be small 
depending on the tunneling time and/or pulse duration. 
However, it is not the case here, as shown in Fig.|21 where 
the different Eq. [31 contributions to the photocurrent are 
shown. Blue squares show the results without the contri- 
bution during the pulse (the second integral in Eq.|Sl) and 
neglecting the biexciton state in the entire simulation. 
The yellow circles show the same simulation but includ- 
ing the extra integral term. The difference between the 
two traces is noticeable and cannot be ignored. Including 
the biexciton in the simulation and using the full Eq. [SI 
results in a photocurrent trace that exhibits decay of the 
oscillations, but increases slightly with the pulse area, 
which is not consistent with experimental observations. 

Experimental results exhibit a decay of the oscillation 
with increasing pulse area, which cannot be obtained 
from our model so far. Notice that any kind of addi- 
tional constant decay or dephasing cannot explain the 
experimental trace either, since increasing the pulse area 
makes this quantity smaller compared with the pulse 
intensity, and the Rabi oscillations are enhanced over- 
all. This behavior indicates the presence of other lev- 
els that contribute to decoherence in the systemiSS It 
has been recently demonstrated that the level contin- 
uum in the wetting layer (WL) plays an important role 
in the background absorption in self-assembled QDs at 
high energies^li We here include the dephasing produced 
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FIG. 3: Photocurrent signal as function of pulse area O for 
1 ps pulse showing different contributions of Eq. Blue 
squares show results neglecting biexcitons and without the 
second integral in Eq.|3 yellow circles include the latter. Tri- 
angles show results including biexcitons and all terms in Eq. 
|5] Dashed line at 7pc = /g ~ 13.1 pA indicates optimum 
value for photocurrent if only exciton channel is considered. 



by the WL continuum, which in essence is being pop- 
ulated non-resonantly by the light pulse and then con- 
tributes to the dephasing of all the low-energy excita- 
tions. It is possible to excite electrons from a bound state 
in the QD (valence band) to the continuum of the WL 
(conduction) and from the continuum (valence) to bound 
(conduction) We can model the WL as a unique level 
detuned from the exciton states by S^j, and coherently 
connected to the ground state by a dipole interaction . 
This results in an additional term to the Hamiltonian Hq 
given by — 6w\w){w\ — ^r2u,(t)|0) {w\ + h.c.. The elec- 
tron or hole in this WL state will scatter quickly and 
leave behind the other localized single particle s. We de- 
scribe this process by one decay rate r^,, which connects 
the WL states to s. This process is shown in Fig. Efa) 
where it is represented by a one-direction arrow. Even- 
tually the electron or hole in this s state will tunnel and 
the system will return to its ground state. This connec- 
tion is modelled by Eg, and given by Eq.[31 Following the 
same analysis as for the biexciton, we would not expect 
to populate these levels for low intensity, since they are 
apparently far detuned from the laser energy. However, 
the WL level broadening allows this process to take place, 
resulting in these extra levels being populated with in- 
creasing pulse intensity, and a growing background signal 
in the photocurrent (as that subtracted in experiments). 
This is clearly shown in Fig.^b), where we compare the 
occupation of some states of the system, with Sw = 20 
meV, = 40 meV, fl^ ^n^, and, F, = F'/" given by 
Eq. [3 Figure 2Jc) shows the resulting photocurrent ob- 
tained from Eq. [SI without the background contribution 
of the WL. The photocurrent exhibits a decay of the oscil- 
lation induced by the nonresonant excitation to the WL 
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FIG. 4: (a) Level structure and couplings including nonreso- 
nant excitation to wetting layer WL. (b) Occupation proba- 
bility of states considered in model as function of pulse area 
Q. (c) Photocurrent signal vs. O for 1 ps pulse. Squares show 
results of our model and circles are experimental data from 
Ref. Q. Notice this model reproduces experimental results. 



and it fits quite well the experimental result (although 
no fine tuning of parameters has been done).™ 

We have described the WL by an additional level that 
is coherently excited out of resonance. In reality, the WL 
has a broad continuous distribution of levels resulting in 
an incoherent pump excitation, since we have a one-way 
excitation, as can be seen in Fig. El[a,b) (the population 



of state s only increases). Therefore, we can propose 
an equivalent model that involves a leakage rate Fwl 
connecting the ground state 

the summation is over all levels that compose the WL. 
This can be written as 



0) to the WL continuum 

2 

S{Eiy — Eq — haj), where 
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where pwL is the WL density of states, and /i^, is the 
effective dipole moment connecting the WL continuum to 
the dot ground state. The master equation ||2Jl with rate 
® gives quantitatively similar results to those shown in 
Fig. He) if we use a density of states pwL — ' 
where 6w is the WL detuning with the laser, and F^ is 
the broadening of the WL levels from our previous model. 
Notice that this channel naturally results in an intensity 
dependent decoherence, as that described recently.^-^^ 

We have presented a realistic model to describe the 
Rabi oscillations observed in a single QD photodiode in 
the presence of pulsed light. Using a model that includes 
multi-exciton and WL states, a density matrix formalism 
shows that the two-level model breaks down and gives 
an imprecise interpretation of the results. Full inclusion 
of the biexciton allowed us to derive an expression for 
the photocurrent signal that corrects the simplest model 
used in the literature, and more precisely assesses the ef- 
ficiency of the level rotation. Our study of the damping 
mechanisms makes us conclude that the only possibility 
to explain the damping with increasing pulse area (and 
fixed pulse duration) is an off-resonant excitation to a dif- 
ferent level. This is identified most likely as the WL con- 
tinuum of levels. This would also explain the background 
in the photocurrent signal observed in experiments pifii 
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